Probing neutrino magnetic moment and unparticle interactions 

with Borexino 



Daniele Montanino 1 , Marco Picariello 1 , Joao Pulido 2 

1 Dipartimento di Fisica, 
Universitd del Salento and Sezione INFN di Lecce 
Via Arnesano, 1-73100 Lecce, Italy 
2 Centro de Fisica Teorica das Particulas, 
Departamento de Fisica 
Instituto Superior Tecnico 1049-001 Lisboa, Portugal 

Abstract 

We discuss the limits on the neutrino magnetic moment and hypothetical interactions with a 
hidden unparticle sector, coming from the first neutrino data release of the Borexino experiment. 
The observed spectrum in Borexino depends weakly on the solar model used in the analysis, since 
most of the signal comes from the mono-energetic 7 Be neutrinos. This fact allows us to calibrate 
the v-e scattering cross section through the spectral shape. In this way, we have derived a limit on 
the magnetic moment for the neutrinos coming from the Sun (in which a v„ and u T component is 
present): \x v < 8.4 x 10 _11 /x^ (90%CL) which is comparable with those obtained from low energy 
reactor experiments. Moreover, we improve the previous upper limit on magnetic moment of the 
v T by three orders of magnitude and the limit on the coupling constant of the neutrino with a 
hidden unparticle sector. 

PACS numbers: 26.65.+t, 14.60.Lm, 13.15.+g, 14.60.St 



1 



I. INTRODUCTION 



The present experimental world neutrino data (apart the LSND experiment which was 
not confirmed by the recent MiniBooNE result) provide a robust interpretation in terms of 
a three active neutrino oscillation scenario (for a recent review, see e.g. [l|). In particular, 
solar neutrinos together with KamLAND data, can be explained in a simplified two-neutrino 
framework (in the limit of a vanishing # 13 ) with 0]: 

Am? 2 = (7.58±8:|S) x l(T 5 eV 2 , 

tan 2 9 12 = 0.56 ±° f 9 . (1) 

However, there is still room for subleading non standard neutrino interactions in the interpre- 
tation of data. For example, an evidence for time modulation of the oscillation probability 
in the solar neutrinos could be explained in terms of a non-zero magnetic moment which 
could induce transitions among active and sterile neutrinos (3)]. 

Recently, the Borexino collaboration has released the first data relative to about three 
months of data taking Q]. In this experiment, solar neutrinos (mainly, those coming from the 
7 Be source) are detected through v-e scattering and the recoil electron energy is measured 
through a scintillation technique. The observed event rate is essentially consistent with the 
one predicted by the Standard Solar Model jl] and the oscillation hypothesis. However, 
since the main source observed in Borexino is the monoenergetic 863 keV 7 Be source, a 
precise calibration of the differential da/dT e (E u ,T e ) cross section (where E u is the incident 
neutrino energy and T e the recoil electron kinetic energy) is possible through a spectral shape 
analysis. For example, neutrino electromagnetic form factors would influence the scattering 
cross section [fjj]. In particular, a non zero neutrino magnetic moment introduces a term 
which grows with the inverse of both the energy of the incident neutrino and with that of 
the recoil electron. For this reason, a low energy experiment (such as Borexino), is in a 
favorable situation. 

Although the limits obtained are still weaker than those obtained by a direct measure- 
ment of the v e -e scattering in reactor experiments [3, [1] (and those that could be obtained 
by Borexino itself from a calibration experiment with an external source of neutrinos or an- 
tineutrinos (9|) we should stress that these are short baseline experiments, and they measure 
the magnetic moment of the v e component. Instead, solar neutrinos embed also a compo- 
nent of and u T , for which the limits are much weaker [l(J ITlf . SuperKamiokande data 
were used in the past to study the neutrino magnetic moment [l2|, [Lj, UM, lis], but since 
this experiment observes the continuous 8 B source, it is difficult to disentangle the effects 
of a spectral distortion due to non-standard interactions and those due to the oscillation 
mechanism. In the past, the Borexino collaboration tried to put a limit on /ij, using the 
prototype of the Borexino detector (the Counting Test Facility, CTF) [13]. However at that 
time, due to the smallness of the CTF, no solar neutrinos were observed and a limit was 
established assuming the theoretical SSM neutrino flux. 

Concerning other non-standard interactions, the possibility of a conformal hidden sector, 
called "unparticle sector" , which couples to the various gauge and matter fields of the SM 
through non-renormalizable interactions has been recently proposed [l8|. The unparticle 
sector is assumed to have a non-trivial infrared fixed point, An, below which the sector has 
a scale invariance and the hidden operators become an effective unparticle operator with 
non-integral scaling dimension d. Limits from low-energy neutrino-electron scattering in the 
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unparticle physics framework have been recently obtained in 19]. We will show that the 
limits coming from Borexino are stronger. 

The plan of the paper is the following. In Sec. [Ill we review the contribution to the cross 
section due to the neutrino magnetic moment or the coupling with the unparticle sector; 
in Sec. Hill we describe briefly the experimental input; in Sec. [IV] we describe our analysis 
technique; in Sec [V] we derive the upper bounds on the neutrino magnetic moment from 
SK and Borexino experiments, emphasizing the fact that the former is quite sensitive to 
the Solar Model assumed while the latter is independent; finally, in Sec. [VI] we draw our 
conclusions. 



II. ELECTRON-NEUTRINO CROSS SECTION 



A. Neutrino magnetic moment 



For a neutrino with flavor a, the scattering process standard differential v a -e cross section 
as a function of the incident neutrino energy E v and of the recoil electron kinetic energy T e , 
is given by [20I]: 1 
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with cr = G 2 F m 2 e / '(2ir) . For /i and r neutrinos, where only neutral current interactions are 
possible, the Standard Model of electroweak interactions provides g v ' T = 2 sin 2 9w — \ and 
g A T = — 1/2, with sin 2 9\y = 0.23122 10]. For electron neutrinos, where also charge current 
interactions are possible, we have g VA — > g v ' T A + 1. 

Besides standard interactions, neutrinos can couple with photons through a possible 
magnetic dipole and/or charge radius. The effective low-energy 1/-7 interaction vertex is 
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where (r 2 ) afc is the charge radius and /i a b and d a b are the neutrino magnetic and electric dipole 
moments respectively. Since for ultrarelativistic neutrinos it is not possible to distinguish 
between the two dipole moments, for simplicity we consider only magnetic moments. For 
an individual incoming neutrino with flavor a, since the outcoming neutrino flavour is in 
general not observable, the only phenomenologically relevant parameter is a combination of 
the magnetic moment matrix /x a &: 



/5> 



ab 



(4) 



For Dirac neutrinos fi a b is a generic complex matrix and involves transitions among left and 
right (sterile) states. Conversely, for Majorana neutrinos the transitions are among neutrino 



1 In this work for completeness we have also included the 1-loop corrections to this formula 2l|. However, 
the effect of such corrections is negligible. 
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and antineutrino states of different flavours. In this case the matrix // {, is antisymmetric 
(and in particular, fi aa = 0). 

Actually, the charge radius could be absorbed in a redefinition of the gy's: 

9v^9v + ^{rl)a^e w . (5) 
However, as we will comment later, Borexino is largely insensitive to variations on the axial 



and vector couplings g v and g\ within the limits coming from present phenomenology [22 
For this reason, we limit our analysis to the magnetic moment assuming standard values for 
g v and g%. 

Since the neutrino flavor composition in the experiments considered here depends on the 
energy range of each experiment, the bounds we derive in this work are actually applicable 
to an effective neutrino magnetic moment, which is a linear combination of the individual 
flavor magnetic moments whose coefficients depends, as will be seen later, on the weighted 
average survival probability for each experiment. 

Robust cosmological arguments show that fi u should be smaller than 10 -8 Bohr mag- 
netons, /is, [111], although other astrophysical arguments largely override this bound (see, 
e.g., [24( and reference therein). However, such arguments are model dependent and thus 
less reliable. The strongest direct bounds on the v e come from the TEXONO experiment 
0], i.e., n e < 0.74 x 10 -10 /i-B and, more recently, from the GEMMA experiment, i.e., 
/i e < 0.58 x 10" 10 /zb at 90% C.L. [8J. However the limits for the /^ jT are much weaker 



(AV < 6.8 x 10- 1( Vb [25j, fj, T < 3900 x 10- 10 /ib [26j) 



The contribution to the v a -e cross section due to the neutrino magnetic moment interac- 
tion is given by [6j: 
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where we have explicitly factorized out the \i a dependence from the expression of the cross 
section. A comment is in order. In principle an interference term between the magnetic 
moment and the weak interaction is possible. However, this interference term vanishes if the 
neutrino are longitudinally polarized and the electrons are unpolarized. If neutrinos cross a 
strong magnetic field (such as the solar one) they can acquire a transverse polarization due 
to the precession. In this case an interference effect could contribute to the cross section 
0. Here we neglect this effect since the first Borexino data was taken in a period of low 
magnetic field activity. 



B. Coupling with Unparticles 

Recently, a scale invariant ( "unparticle" ) sector which decouples at high energy was 
proposed in [l8|. Leptons can couple for example with a scalar unparticle sector 18, 27, 2^ 
through the Lagrangian 2 

£u = Xe-^TeO u e + £ Xf-^u^u^ + h.c. , (7) 

A U a,b A U 



2 In general it is possible also the case in which a coupling with a vector or tensor unparticle is considered. 
Here, for simplicity, we consider only the scalar case. 
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where d is the non-integral scaling mass dimension of the unparticle operator, Ay is a typical 
scale of the unparticles physics and it can be assumed ~ O(TeV), Ou is the unparticle 
operator, and the A's are the coupling constants of the leptons to the unparticle sector 
(possible flavor changing interactions v a — > have been also taken into account). The 
contribution to the scattering amplitude for elastic u a -e scattering from the exchange with 



a scalar unparticle is [19 



M* = [Hkf>a{h)\ j— ^ [e(p f )e( Pl )} , (8) 

where q = kf — ki and 



Itt 5 ' 2 r(d + l/2) 
(27r) M sin(7rd)r(d-l)r(2d) 

From this amplitude, the contribution to the v a -e cross section can be calculated: 3 



K-^r(K, T e ) = X 2 a vE2A «U {m e T e ) M -\T e + 2m e ) , (10) 



and we have conveniently defined 



^=^/E(IWe) 2 . (11) 

As in the case of magnetic moment, the final state will be a sterile right state or an antineu- 
trino state depending on the Dirac or Majorana nature of the neutrinos. Notice that for 
d < 3/2 the cross section diverges for low T e , thus low energy experiments are most sensitive 
to the unparticles. Notice that for d — 1 we have the same T~ l dependence as in the case 
of the magnetic moment. 



III. THE EXPERIMENTAL INPUT 

The Borexino experiment at the Gran Sasso underground laboratory is designed to 
study mainly the 863 keV monoenergetic 7 Be solar neutrinos, through a real-time and low- 
background detector. A detailed description of the experimental apparatus and its ancillary 
plants can be found in [29T ]. Briefly, Borexino consists of 300 tons of an high purity liquid 
organic scintillator (pseudocumene, C9H12, p=0.88 g/cm 3 ) doped with PPO at 1.5 g/1. The 
scintillator mixture is contained inside a thin nylon sphere (8.5 m in diameter). This vol- 
ume is viewed by 2212 8" photomultipliers (PMTs) installed on a stainless steel sphere of 
13.7 m in diameter. The scintillator is shielded against background from PMTs and other 
external sources by a 2.5 m buffer of pseudocumene, an outer nylon vessel and 2 m of high 
purity water. For solar neutrinos search a fiducial volume of 100 tons is selected off-line. 
In Borexino solar neutrinos are detected through the scattering reaction v + e — > v + e. 
The recoil electron energy is converted into light inside the scintillator. The intrinsic 14 C 



3 One can easily check that for almost massless neutrinos the interference term with the standard amplitude 
is negligible. 
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contamination 4 and finite energy resolution set the detection threshold at about 200 keV. In 
Borexino no directionality is possible to search for neutrino interactions and it is not possible 
to distinguish on an event-by-event basis between neutrino processes and fi/j backgrounds. 
Therefore, the radiopurity of the scintillator is a fundamental experimental issue. The first 
results 0] have shown that the radiopurity achieved is beyond the expectations and this 
allows to extend the research program as we attempt to do in this work. 

The data taking started in May 2007. The collaboration released the first data in Au- 
gust 2007 The observed flux of 7 Be neutrinos is (within the errors) consistent with 
the Standard Solar Model prediction in the hypothesis of Mikheyeev-Smirnov-Wolfenstain 
oscillations, i.e., 47 ± 7(stat.)±12(sys.) counts per day (cpd) in 100 tons for the 7 Be 
(863 keV) neutrinos, against a theoretical value of 49±4 cpd's. The collaboration has also 
produced a spectrum of the observed events with the visible electron energy K e in the range 
270 < K e < 800 keV. This spectrum is shown in Fig. 6 of [H, in which the number of 
events per day and per 100 tons of scintillator are shown in 53 bins of energy (we show our 
equivalent plot in Fig. [T]). The collaboration quote a 15% uncorrelated error for each bin. 
However, since there are also bins with less of three events, the mere statistical error for 
these bins would be greater than 15%. For this reason we prefer to be conservative and sum 

— 1/2 

in quadrature the 15% error quoted by the collaboration to the statistical error N { , where 
Ni is the number of events in each bin. In the majority of bins, where N{ 3> 3, the statistical 
error is negligible and the uncorrelated error is just that quoted by the collaboration. The 
main source of correlated error comes from the determination of the fiducial mass. The 
collaboration quote a 25% error equally correlated among all bins. 

From |4| we know that the measured light yield is about 500 photoelectrons/MeV. The 
light yield affects the energy resolution of the detector and at a first approximation we can 



assume a gaussian energy smearing with a = y T e /500. We notice that this assumption does 
not work well in the low energy regime (mainly below 200 keV) where non linear effects due 
to quenching take place. In the energy range we are considering, namely [270,800] keV this 
effect is expected to be on the order of a few %'s. With the above assumption the resolution 
function for the detection of solar neutrinos is thus: 



K(K e ,T e ) = exp 
v lira 



(K e - T e f 
2a 2 



(12) 



where T e is the real kinetic electron energy, K e is the visible (measured) one, and o = 

The main sources of background in the detector come from the decay of contaminants 
contained in the scintillator (being the cosmogenic contribution almost completely rejected 
by muon vetoing and other techniques with the exception of n C which, however, could 
dominate the spectrum above 1 MeV). Internal background is due to the (3 decays of 210 Bi 
and 85 Kr and the a decay of 210 Po as reported in [i| . The former background is enormously 
reduced by a//3 pulse shape discrimination (PSD). However, the rejection efficiency of the 
PSD is not 100% and may change with the energy. In our analysis we adopt the same 
assumption used in j3] and take into account a possible small residual of a-like events in the 
/3-like/neutrino events spectrum by means of a gaussian peaked around 410 KeV: 

S a (K e ) = A^iopo • K(K e , 410KeV) . (13) 



The C have a f3 decay with an end-point energy of 156 keV. 



6 



Here A/2iop is the unknown normalization to the spectrum that should be determined by 
the fit. 

As pointed out above (3 decays cannot be rejected. They contribute to the total spectrum 
through the 

S P {T e ) = E^(l-|) 2 || F B {T e ) , (14) 



where B E { 210 Bi, 85 Kr}, E e = T e + m e and p e = ^J2m e T e + T 2 are the total electron energy 
and momentum, and the Fermi correction Fb is given by [30]: 



F B (T e ) = 2(1 + 7o ) (2R Pe )^- 2 e- , (15) 

where R = 0.426y4 1//3 a e m Jm e is approximatively the radius of the nucleus and, in turn, 

v(T e ,Z) = a e . m ZE e /p e , 



7o = ^l-a 2 m Z 2 , 

with A and Z atomic number and charge of the parent nucleus respectively. The Q values for 
the two f3 sources are 1162.1 KeV for the 210 Bi and 687.1 KeV for the 85 Kr. Moreover, since 
the (3 decay of the 85 Kr is forbidden, the spectrum should be corrected by the multiplicative 
factor pl/Q 2 + (1 — T e /Q 2 ) 2 [3(|. Of course the observed (3 spectrum are the convolution of 
the true spectrum with the resolution function. 



IV. ANALYSIS 

The observed total spectrum in Borexino is given by 

S(K e ) = EM<D S / dE u ip s (E u ) J2 K e ) ■ P ea (E u ) + S a (K e ) + S p {K e ) , (16) 

s J a die 

where s G {pp, pep, 7 Be, 16 0, 14 N} are the solar sources (the contribution from 8 B and 
hep neutrinos is negligible) with their (normalized) spectra tp s and Standard Solar Model 
flux $ s (in particular, we have used the AGS05 fluxes, see Table 6 of [5]), P ea (E u ) is the 
oscillation probability P{y e — * u a ), and J\f s are normalization factors (with the constraint 
J\f s > 0); the "tilde" means that the cross sections/background spectra have been convoluted 
with the resolution function. In practice the CNO sources (i.e., the 16 0, 14 N) give a small 
contribution to the spectrum. Moreover, their spectra are almost indistinguishable among 
them and from those of the 210 Bi background. Since we do not expect a strong spectral 
distortion from oscillations, following the Borexino paper, we have embedded the CNO 
contribution into the 210 Bi spectrum. 

For 7 Be and pep neutrinos the spectrum is almost monoenergetic apart a few keV broad- 



ening due mainly to collision and thermal effects [3jj. In the energy range explored by the 
experiment the only contribution to the spectrum comes from the 863 keV 7 Be line, the 
contribution from the other line (385 KeV) being well below the threshold. In our analysis 
we fix only the pp and pep neutrinos at their Standard Solar Model value (i.e., we force 
A/"p Pi p e p = 1) [5], as their sources are affected by a very low theoretical uncertainty (~l-2%). 
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All the other solar fluxes, the background normalizations and, of course, the non-standard 
parameters, are taken as free variables in the fit. 

Regarding the oscillation probability, we have used the standard Mikheyev-Smirnov- 
Wolfenstein (MSW) probability with the oscillation parameters defined in Eq. (CQ). 5 We 
do not fit the mass-mixing parameters with Borexino data since it is not the goal of this 
work. A comment is in order. In principle when non standard interactions are present, 
also the oscillation probability could be affected by new interactions. In particular, neutrino 



magnetic moment can induce resonant spin-flip precession 33| (for a recent review, see also 



34(, Sec. 13.1 and reference therein). Anyway, barring oscillations into sterile neutrinos, 
any change in the probability would reflect only in a very small change of the spectral shape 
(since it is determined mainly by the 7 Be line), while the overall normalization is left free. 6 
For this reason the precise choice of the probability is not critical to our scope, provided 
that: 1) no active-sterile oscillations are allowed and 2) no dramatic spectral distortion in 
the CNO spectra are expected. 

We now have all the ingredients to calculate the x 2 '- 7 

X 2 = E(^ h - N? xp )(°- 2 UN^ - Af xp ) , (17) 

y 

where N^ xp (Af Th ) is the experimental (theoretical) number of events in the i-th bin, and, 
as explained in the previous section, the matrix cr?- is calculated as: 

a 2 = [Af xp + (0.15Af xp )% + (0.25iVf xp ) ■ (0.25Af xp ) . (18) 

Since all the unknown parameters appear linearly in the calculation of N? h , the \ 2 mini- 
mization is straightforward. Defining 

f)N Th 

F " = -k (19) 

(with k G {pp, pep, 7 Be, CNO+ 210 Bi, 85 Kr, 210 Po}), it is easy to show that 

A4 = E [(Fa- 2 ^)- 1 ^- 2 ] ki Uf xp - £ M s fA , (20) 

i \ se{pp,pcp} / 

where the pp and pep sources are fixed and not fitted. In any case, negative A4's are not 
allowed. 

The "best fit" spectrum with standard interactions only is shown in Fig. [TJ which is 
similar to Fig. 6 of [4|, with black solid line. (In the figure is shown also the spectrum with 
a non zero magnetic moment with the black dashed line). For comparison, we also report 



We have used the approximate formulae for calculating the P ee survival probability already averaged over 



the production zone given in 32 1. 

6 Moreover, to simplify our analysis, we do not allow any v «-> v transitions so that z7's are absent in the 
flux. 

7 In presence of bins with zero or few events the x 2 function should in principle be corrected as prescribed in 
jlOj. However, since the low statistic bins are less relevant for the analysis, we prefer to use the standard 
X 2 ■ In this way the unknown parameters can be extracted analytically. 
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FIG. 1: (color online) The "best fit" spectrum with zero (black solid line) and nonzero (black 
dashed line) magnetic moment. The contributions of the 7 Be source (red medium-dashed line) 
and those from the CNO+ 210 Bi ( green dotted-dashed line), 85 Kr (blue dotted line), and 210 Po 
(light blue log-dashed line) in the case \x v = are also shown. The number of counts per day in 
the "best fit" (fi v = 0) case in the whole energy range are 49 for the 7 Be, 12 for the CNO+ 210 Bi, 
18 for the 85 Kr, and 1 for the 210 Po source. For illustration, we also show the contribution coming 
from the magnetic moment only for /i„ = 7x 10 -11 //£ (red dotted line). 



the number of events per day in the full recoil energy range for each source for 100 Tons of 
scintillator in our "best fit" case: 49 for the 7 Be, 12 for the CNO+ 210 Bi, 18 for the 85 Kr, 
and 1 for the 210 Po. 

The value of Xmin * s 37.6, is slightly lower than the one obtained by the collaboration 
(Xmin = 41.9), due to different assumptions in the two analyses. In the same figure the 
various contributions (apart those from from pp and pep neutrinos which are very small) 
are also shown. The main contributions to the spectrum come from the CNO+ 210 Bi, 85 Kr, 
and 210 Po. Among these, only the first could be slightly affected by the functional form of 
P ea (E u ). We see also that the this contribution as almost flat, while those from 210 Po is a 
peculiar "bump". The main contribution for the spectral distorsion at low energies (thus 
mimicking those coming from non-standard interactions) comes from the 85 Kr background. 
As we will discuss later, this background is the main limitation to the measure. 



V. LIMITS ON THE NON STANDARD INTERACTIONS 



Introducing the non standard interactions, we see that the main contribution to the 
spectral distortion comes from 7 Be neutrinos. In particular, we see that, if only 7 Be neutrinos 
are taken into account, the contribution to the spectrum coming from the non standard 
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interactions is: 

J~ non std 

5S(K e ) = &.-? IR -(E v ,T e ), (21) 

where £ = /i u ,X u and da Don > std /dK e is the non standard contribution given by Eqs. (El) or 
ffTUj) . after the proper convolution with the resolution function. The label "eff" means that 
we must consider an effective coupling, given by: 

& = £ p «(£S)-£, (22) 

a 

where E® = 863 keV is the energy of the 7 Be neutrinos. Due to the unitarity of the prob- 
ability, this is also equivalent to introducing an equal magnetic moment (/x„) or unparticle 
coupling (Aj,) to all flavors. Since also non-monochromatic sources contribute to the spec- 
tra, for simplicity we stick into this simplified hypothesis of equal non-standard parameters. 
However, since the final result of the analysis does not depend critically on the functional 
form of P ea , the limits on \i v and X u are practically a limit on // e g and X c s. 

In particular, if we trust that the oscillation is simply given by the MSW effect with small 
6*i3, we have that /i 2 ff and A 2 ff are given by 

/4r = P 2v {K) ■/£ + [!- P 2v {ED] (cos 2 6 23 ■ + sin 2 9 23 ■ ^) (23) 

(and, of course, a similar expression holds for A e ff) ■ In principle, in the previous equation 
also interference terms among oscillation amplitudes can be present. However, with the 
oscillation parameter given in Eq. (TjQ) these terms average out and can be safely neglected 
@, EE). 8 In Eq. P 2v {El) is the two-neutrino P{u e -> u e ) standard MSW probability 
and sin 2 9o 3 = 0.38 -~ 0.63 (2a) from recent atmospheric and accelerator (K2K and MINOS) 
analysis pj. Moreover, for low energy neutrinos we have also that with good approximation 
P 2u ~ cos 2 #12 ~ 1/2. 



A. Limits on the magnetic moment 

We start by deriving the upper bound on the effective magnetic moment from the Su- 
perKamiokande (SK) experiment. To this end, the flux measured by SK is: 



/ dE v (ps B {E v ) Ea Pea(E v ) af d {E u ) + tia»(E v ) 

$SK = 



/ dE„ips B (E u ) Ea P ea (E u )af d (E u 



(24) 



where <r^ td and a 11 are the total (standard and non-standard) cross sections in the SK energy 
range. In this case is not possible to define an effective magnetic moment as in Eq. (1231) . so, 
for simplicity, we have assumed an equal \i v for all flavors. 

From Eq. (f24j) we realize that the bound on \x v depends on the total 8 B flux (and hence is 
solar model dependent). For example, using the AGS05 model where $sb = 4.51 x 10 6 x 
(1±0.12) cm-V 1 with the SK measured rate R SK = [2.35 ±0.02 (stat.)±0.08 (sysj] x 10 6 
cm- 2 s" 1 we obtain ^ < 2.1 x IO^Vb (90% CL). On the other hand for the GS98 model 
with the higher 8 B flux $s B = 5.69 x 10 6 x (1 ± 0.16) cm" 2 s _1 we obtain a stronger bound, 



Note that our definition of /i a is different from those used in ref. 15 1 
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FIG. 2: The x 2 as a function of the neutrino magnetic moment /i for standard axial and vectorial 
coupling. We also report the 90% C.L. (x 2 - Xmin = 2.71). 



\i v < 1.3 x 10~ 10 fi B (90% CL), which is more in agreement with the bound obtained by the 
spectral analysis done by the SK collaboration [l6j |. 

As far as Borexino is concerned, the corresponding upper bound on the neutrino mag- 
netic moment is not dependent on the Standard Solar Model assumed. In fact the 7 Be 
normalization is extracted from the experimental data. In order to show the effect of the 
magnetic moment, in Fig. [I] in black dashed line we plot also the theoretical spectrum for 
\i v = 10~ 10 /ib (which is well beyond the 90% limit). From the plot we see that (as ex- 
pected) the spectrum grows at low energies. This behaviour, as we discuss later, can be also 
mimicked by the 85 Kr background. 

In Fig. [2] we show the \ 2 as function of the neutrino magnetic moment where all the 
free fluxes and backgrounds have been marginalized. We see that the 90% C.L. limit on 
the magnetic moment is \i v < 8.4 x 10~ n /iB- 9 We have also tried to fit the neutrino charge 
radius and more in general the neutrino vector and axial couplings gy and gA, as proposed in 
35| . Unfortunately the limits obtained with Borexino are far from being competitive from 
those obtained by other experiments [22j. However, we have verified that, varying gy and 
gA inside the allowed region(s) in Fig. 2 of [22] our limit on fi u does not vary appreciably. 

Although this limit is less competitive than those obtained from reactor experiments like 
in GEMMA (/i e < 5.8 x 10~ n /iB at 90% C.L., Q), we should bear in mind that reactor 
experiments are short baseline. For this reason, the limits obtained in these experiments are 
essentially bounds on the \i e component. Instead, the limits coming from solar neutrinos 



9 In |15| a perspective analysis of the Borexino experiment was performed. Although they obtained a more 
stringent limit on a combination of the Majorana magnetic transition moments, they assumed a fixed 
background. As we will see later, the main source of uncertainty comes from the lack of knowledge of the 
true number of background events. 
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can be translated into limits for jj^ and fi T . In particular, we have seen that for Borexino 
the limit on \x v is practically a limit on fi c g. Using Eq. (|23|) and P 2u ~ 1/2 we obtain a 
conservative limit on \x T (obtained in the worst = and sin 2 612, cos 2 #23 taken 

at their 2a minimum allowed values): 

H T < 1.9 x 1(T 10 /ib, (25) 

which is three order of magnitude stronger than those quoted by the Particle Data Group 



(/i T < 3900 x 10 10 hb jlOQ). Equivalently, we get the limit 

/V < 1.5 x 10- 10 fi B . (26) 

The "plateau" in the \ 2 f° r ~ 0.85 x 10~ 10 /iB is due to the partial compensation 
between the 85 Kr background and the magnetic moment contribution (we have also checked 
that the other normalization factors are almost insensitive to the value of fj, v ). In fact both 
the magnetic moment and the 85 Kr, have the same effect in the spectrum, i.e., to increase the 
slope at low energies. Increasing \i v the slope is kept almost constant if one simultaneously 
decreases the normalization of the 85 Kr background. For example, this can be also seen in 
Fig. [1] from the comparison between the red dotted line (which is the contribution of the 
magnetic moment to the electron spectrum due to 7 Be neutrinos, for \i v = 7 x 10^ u /i B ) 
and the blue dotted line (which is the contribution due to the 85 Kr in absence of magnetic 
moment). At \i v ~ 0.85 x 10 _10 /zb the 85 Kr normalization vanishes and compensation is no 
longer possible. This can also be noticed in the abrupt discontinuity of the derivative of the 
X 2 - 

The main limitation for the measure of //„ thus comes, not from the limited statistics, but 
from the imprecise knowledge of the 85 Kr background. If in the future this background will 
be reduced by a purification campaign, 10 the limit on \x v would become much stronger. We 
also mention that by improving a lot the exposure it might be possible to measure the 85 Kr 
contamination by means of correlated events as reported in [4]]. Moreover, due to the Earth 
orbital eccentricity, the solar flux has periodical variations in one year, while the internal 
background is expected to be nearly constant in time. This means that a better limit on \x v 
could be obtained simply measuring the spectrum in different periods of the year. 



B. Limits on unparticle coupling 

We have done the analysis done in the previous section but using the unparticle cross 
section (flOl) . All the considerations made in the previous section (in particular, the partial 
compensation between the 85 Kr background and the unparticle contribution) apply also in 
this case. 

In Table |T] we show our limits on \ v as function of the dimension d. We fix the value of 
Ay to 1 TeV. Of course, if we choose a different value for Ay, the limits on should be 
rescaled according to Eq. (fit)]) . The values of d (first column) have been chosen in order to 
make a comparison with the limit obtained in Table 1 of [19]. In the third column we show 
our 90% C.L. limits on X u obtained in the same way of the previous section, i.e., taking 



We notice that due to the long mean-life of the 85 Kr, the fate and impact of this background is different 
than that of 210 Po. 
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10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 1 10 10 2 

A y /(Au/TeV) 2d " 2 



FIG. 3: The 90% C.L. allowed zone in the plane (d, A„) (for A u = 1 TeV) as from Ref. (light 
gray), with unconstrained 7 Be flux (dark gray) and with the 7 Be flux fixed by the Standard Solar 
Model. See the text for more details. 



the 7 Be solar flux completely unconstrained. In the fourth column we show the same limit, 
but taking into account the 10.5% theoretical uncertainty on the 7 Be flux This is done 
simply adding a penalty function to the \ 2 

xW+ (^Lziy. (27) 



V 0"7 Be 

As expected, in this case the bounds on X u slightly improve. 11 In Fig. [3] we show also the 
90% C.L. allowed zone in the plane (d, X u ) as from Ref. [l9|] (light gray), with unconstrained 
(dark gray) and constrained (black) 7 Be flux. 

We see that our limits are even stronger than those obtained by Balantekin and Ozansoy 
(l9| . We remark that the limits obtained in [l9j are obtained using short baseline reactor 
data thus sensitive only to z/ e , while our result applies to the combination A pff i n Eq. (1221) . 
Since P ee ~ 1/2, sometimes our limits on A e alone can be weaker than those in r|. However, 



combining our limits on A e g- with those on A e in [19j we can obtain for the first time bounds 
on the coupling constants A M and A r according with Eq. 



We have checked that in the case of the magnetic moment the improvement is negligible. 
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d K (Table 1 of [19]) \ v (unconst. 7 Be) \ v (SSM 7 Be) 

1.01 3.5 x 10" 6 2.4 x 10" 6 1.8 x 1(T 6 

1.05 7.3 x 10~ 6 4.8 x 10~ 6 3.4 x 10~ 6 

1.1 1.9 x 10~ 5 1.1 x 10" 5 0.8 x 10" 5 

1.2 1.2 x 10" 4 0.6 x 10" 4 0.4 x 10~ 4 

1.3 7.2 x 10" 4 4 x 10" 4 2.9 x 10" 4 

1.4 4.5 x 10~ 3 1.8 x 10" 3 1 x 10 -3 

1.5 2.7 x 10~ 2 1 x 10~ 2 0.5 x 10~ 2 
1.7 9.5 x lO" 1 2.4 x 10" 1 1 x 10" 1 
1.9 24.5 4 2 

TABLE I: Confrontation between our 90% C.L. limits on X u with 7 Be unconstrained (third column) 
and fixed by the SSM fourth column for different value of d, with = 1 TeV. Also the limits 
obtained in 19( are shown (second column). 

A comment is in order. In principle, unparticle interactions could affect also the produc- 
tion and propagation of neutrinos 3f|. However, we have already stressed that the major 
contribution of the solar flux comes from the 7 Be monoenergetic line. This largely compen- 
sate eventual uncertainties in the conversion probability. For this reason, we think that the 
analysis with the completely unconstrained 7 Be should be considered to be more reliable. 

Recently, in [37] it has also pointed out that a generic unparticle scenario will generate 
contact interactions between the particle operators. This contact term would generate a 
scattering amplitude similar to those in Eq. (jSj), but with d = 2 and T 2 = l. 12 Assuming 
that this term is dominant (and thus any interference term is negligible), the v-e cross section 
is thus similar to those in Eq. (fTUj) but with 



A 



cont 



1/2 
ab\2 



'ev J 

b 



(28) 



where A"* are the contact term couplings for the v a e — > v^e scattering. The equivalent limit 
on A™ nt for Ay = 1 TeV is A cont < 3.5 (1.4) in the hypothesis that 7 Be neutrino flux is 
unconstrained (constrained). We do not perform a combined analysis with \ v and A™ nt 
unconstrained. However, since we see that the limit on is A£° nt in generally much weaker 
than those obtained for d < 2, if the two coupling constants were of the same magnitude, 
the contribution to the cross section due to the non-contact term would be dominant. In 
this hypothesis, the limits on X u in Fig. [3] can be safely assumed. 



VI. CONCLUSIONS 



In this paper we have analyzed the first Borexino data release to constrain the neutrino 
magnetic moment and the coupling of neutrino and electrons with an hypothetical unparticle 
sector. The analysis is performed analyzing the spectrum of the recoil electron energy. 
Since the leading contribution to this spectrum comes from the monoenergetic solar 7 Be 



We thank B. Grinstein for pointing us this fact. 
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neutrinos, the shape of the spectrum is almost independent from the energy dependence 
of the oscillation probability. The other contribution to the spectral shape is due to the 
internal background of the detector. 

We have performed a x 2 At assuming unconstrained both the solar fluxes (except for the 
pp and pep, which, however, give a negligible contribution) and the internal backgrounds. 
In absence of non-standard effects our results are in good agreement with those obtained 
by the Borexino collaboration. 

An upper limit on the neutrino magnetic moment is found: \x v < 0.85 x lCT 10 /!£. Al- 
though this is not the strongest limit in the literature, we stress that the magnetic moment 
measured by Borexino is a linear combination of the magnetic moments of the different 
neutrino flavors. Within a reasonable assumption on the oscillation probability this limit 
translates into an upper one for the (much more unconstrained) fi T . 

In the same way, we have obtained limits on the couplings with the unparticle sector. The 
results are summarized in Table [H In this case the limits are stronger than those obtained 
in the literature Furthermore, in this case the coupling constant is actually a linear 

combination of the couplings of neutrinos with different flavors. In this way we have obtained 
upper bounds on the still unknown parameters A M and A r . We want here to stress the fact 
that this limit strictly applies if only scalar and contact unparticle operators are present. If 
other operators are included in the analysis the bounds would be relaxed. A more complete 
analysis with further unparticle operators is deserved for a time being. 

Finally, we are confident that in future Borexino will improve these limits. Since the effect 
of non-standard interactions is partly hidden by the internal contaminants, a reduction or 
at least a better knowledge of these backgrounds would be of a great help. With at least one 
year of data-taking, using the seasonal variations of the solar neutrino flux due to the orbital 
eccentricity (about 3% in one year) Borexino will be able to disentangle the contribution to 
the spectrum coming from the 7 Be neutrinos to that due to the (if constant) background. 
Moreover, a substantial reduction of systematics and a dramatic increase in the statistics is 
also expected. 
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